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ABSTRACT 


The earth's atmosphere is considered as made up of oblate spheroidal 
layers of variable density lying over an oblate spheroidal Earth. The 
gravitational attraction of the atmosphere at exterior points is com- 
puted and its contribution to the usual spherical harmonic gravita- 
tional expansion is assessed. The potential is also found for points 
at the botttom of the model atmosphere. This latter result is of 
interest for determination of the potential at the surface of the 
geoid. The atmospheric correction to the geoid determined from ‘ 
satellite coefficients is given. 
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A SIMPLIFIED MODEL FOR THE GRAVITATIONAL 
POTENTIAL OF THE ATMOSPHERE AND ITS EFFECT ON THE GEOID 

by 

Stephen ,j, Madden, Jr. 

I. Introduction 

The gravitational contribution for the earth's atmosphere is be- 
coming of more interest as the precision of gravitational determinations 
increases and as it becomes important to relate satellite altimetry 
measurements to the gravitational field of earth. Accordingly, it is 
useful to have estimates of the importance of its effect. The present 
report is concerned with a simplified treatment which is analytically 
tractable and which, at the same time, preserves the essential features 
of the actual situation. 

Since there is a flattening effect due to the ellipticity of the 
earth, the atmosphere is considered to be made up of oblate spheroidal 
layers of constant density. The layers are confocal 'with the reference 
ellipsoid defined by the adopted reference constants of the 1968 I.A.U. 
The lower limit of the atmosphere is assumed to be the reference 
ellipsoid. 

The general procedure used is to find the contribution, to the 
external and internal potentials, of a thin spheroidal shell of constant 
density. This shell is bounded by confocal ellipsoids. The external 
ellipsoid has semi-major axis a', semi-minor axis c* and eccentricity 
e’; the corresponding quantities for the internal ellipsoid are a, c, 
and e. Since the ellipsoids are confocal, 



Once the potential of the shell is found, a limiting argument is used 



to find the contribution of a shell of finite thickness and variable 
density. 

2 . The E xternal Potential 

The solution to the external problem follows easily from a result 
in (Macmillan, 1958). !£ the field point is located by spherical co- 

ordinates r and 6 outside of a solid ellipsoid of density p. then the 

potential is 


WG a J 0 TKTTtBhW P 2n (coS ">■ 


( 2 ) 


where M is the mass of the spheroid, G is the gravitational constant, 
r is the distance from the center, and 6 is the colatitude. If the 
potential due. to a spheroid of density p and semi-major axis a' is 
computed using (2), and if the potential due to a spheroid of semi- 
major axis a and density -P is computed, then the potential caused by 
the spheroidal shell contained between the two bounding surfaces is 
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n=0 

Here AM is the mass of' the shell, and r > a. The total potential 
contribution of a sequence of very thin layers is, using a limiting 

argument, 

(-D n 
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The external potential of the atmosphere depends only on the total mass 
and the- linear eccentricity, B. This result is in agreement with the 
more' general MacLaurin’s theorem (MacMillan, 1958). In this expression, 
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the quantity E is the same as that for the reference ellipsoid. It can 
be seen, from a comparison test using l P 2n^ c0S < t ^ a ' t this series 
converges at all points on the ellipsoid. 


3. The Internal Potential 

The derivation of the internal potential is more complicated. 

Since there is an axial symmetry in the problem we need only find the 
potential on the axis of symmetry, the z-axis. (MacMillan, 1958, pg. 
360). The first step is to find an expression analogous to (2) where 
the field point lies inside a homogeneous ellipsoid of density p. This 
expression follows as a special case of an expression in (MacMillan, 
1958), for this potential, 


V = 7tpGa 2 c ^ |l 


z 2 ^ ds 

c2+sl V(a 2 + s 3 2 (c 2 + s) 


L 


The integral here can be evaluated' directly in terms of elementary 

✓ 

functions (Pierce and Foster, 1961), or found as a limiting case of 
the triaxial potential (Byrd and Friedman, 1971). It is 
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Here z lies within the ellipsoid. 

If we consider the contribution to the -potential of a thin shell 
between ellipsoids with a’ = a + A and a, where A/a << 1, then if z 
lies inside both shells, 
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If we look at the case where A/a << 1, then 


* 2 - a 2 = 


a = 2aA . 


Similarly if we define 


2 - 1 

f (a) » a c sin e, 


with 


2 f 2 • A „ E 
a -E , e - - , 


then 


f(a’) = f(a) + f»(a)A 


If we carry out the differentiation, 


2 c'sin‘ 1 e’ -a 2 c sin’ 1 e z (2ac sin -1 e + ~~ sin 1 e - Ea) A 


Thus, on the 2 -axis, 


AV 


shell 


- j ( 1 + j [ 2ac sin‘ 1 e + ~~ sin" 1 e-Ea] A- 


2 z‘ 


4 


or 


AV shell = 2,tGpA II + z2 p 


2z 2 a 


with 


1 a 3 -1 

M * 2ac sin e — sin e - Ea. 


( 5 ) 


The total contribution of many thin shells to the potential is more 
cumbersome than in the exterior case since the density variation must 
be explicitly considered. Let us rewrite (5) as 


AV shell ' 2TTGQ(a)p(a)da, 


( 6 ) 


with 


da = A 


and 


Q(a) 


M 2 M _ 2z a 
E E 3 E 2 


If we sura the total contributions of many shells, using (6), 


V ATM - 2 * G Ac«P(Od 5 


where a now corresponds to the reference ellipsoid. For present 
purposes the density will be assumed to be an exponential function 
of the quantity £, 

A;.’ e-a ... 

PiS) = p^e H j 5. 



with H a scale height characteristic of the lower atmosphere where the 
bulk of the atmospheric mass is located. A typical value for H is BKm 

With this assumption, 

_£- a 

V ATM “ 2¥Gp 0 /* Q(C)dC ‘ 

\ 

In the evaluation of this integral we can take advantage of the rela- 
tive sizes of a and H; the quantity QCO varies slowly during a scale 
height variation in £ and, as can be shown, 

V ATM " 2^P 0 Q(a)^ e H «, 


or 

V ATM = 2TlG Po HQCa)> 

where terms on the order of H/a have been neglected with respect to one 
If we return to the notation of (5) , 
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( 7 ) 


This expression can be evaluated explicitly if the reference values 
of the constants are used. However, since we are primarily interested 
in representative values, we can utilize the fact that the eccentricity 
of the reference ellipsoid is small to find that 


M * 2a 2 e + a 2 e 5 + 0(e 7 ) . 
If this is used in (7) , 
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V ATO =2l,p 0 GH [ 2a + T5 z2 H’ 

This expression yields the potential anywhere inside the atmospheric 
shell simply by replacing z 2 with z 2 - (x 2 +y 2 ) /2 , and finally . 
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This equation, plus that of (3) yields the desired potential expressions. 


4. Examination of the Solutions 

The most natural questions concerning the external gravitational 
potential of the atmosphere are those involving its contribution to 
the already measured coefficients of the total external potential of 
the earth and atmosphere. The models of- the atmosphere used here 
possess symmetry with respect to the rotation axis of the earth and 
hence can contribute only zonal terms. We thus restrict consideration 
to these terms. The standard zonal expression is given as 


v T 4|-j 2 Gfw- 6) ! ’ 


or 
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with the subscript T referring to the contribution of both the solid 
earth and atmosphere. If we write the corresponding external contri 
bution to the atmosphere, using (3), 
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With this form we can compare the magnitudes of J 2 and and the at- 
mospheric terms. If we use the subscript A to denote quantities per- 
taining to the atmosphere alone, then 

_ ei M A • 

J 2A " 5 M x 

and 

_ 3e 4 M A 

J 4A ' 3$ M t ■ 

The constants appearing here are the eccentricity, available from the 
I.A.U. flattening, f = 1/298.25, and from (Verniani, 1966) are 

e ° .0066945, 

M a = 5.136xl0 21 gm , 

M t = 5 . 976 x 10 22 gm 

With these values we find 

J 2A = -l.lSlxlO- 9 . _ ' 

J,, = 3.302X10' 11 . 

4A 

This value for J 2A compares favorably with the value -1.9 10 9 arrived 
at by Lundquist (private communication) who assumed all of the atmos- 
pheric mass was concentrated in a shell of the dimensions of the reference 
ellipsoid. These numbers correspond to the total values- 


8 



J 2 = 1082.7x10 
J 4 - -1.649xl0' 6 . 

The atmospheric values enter in the seventh place in Jj and in the 
sixth place in J 4 . 

It is probable that the actual, non- idealized, atmosphere ac- 
tually has a larger influence than those shown here. The atmosphere 
is not of constant density on spheroidal shells, and an examination 
of data shows that the earth probably has a somewhat larger value of 
J 2A than is shown here. 

The case of the interior potential is somewhat different , here 
there is no convenient series for comparison. For this reason we com- 
pute, as a measure of importance, the maximum interior attraction of 
the atmosphere at a point on tho surface of the reference ellipsoid. 
From equation (8) it follows that, at any point inside the atmospheric 
shell, the force- is 
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Now on the reference ellipsoid 


x 2 +y 2 _ + ll = i 

3 2 c 2 ’ 
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so that on this surface 


iw 2 = (Is *»o gh l^)! a2+ ( 4 - 7) 


2 

This last quantity varies as a function of z which has 
2 2 

value z = c . Thus 


I -ATM I 15 7 Tp o GHe a 

max 


If we use the numerical values 


p n = 1.23xl0' 3 gm/cm 3 


G = 6 . 67° x 10" 8 cm 3 /gm sec 2 


H = 8x10° cm 


e* = . 0066945 


£ = i.f - 1 - — 

a 1 298.25 * 


then 


|F atm | = 1.5X10" 6 cm/sec 2 = 1.5 xio -3 mgal 


max 


This result is extremely small. It corresponds to 
displacement of the geoid. ' “ 


the maximum 


a negligible 
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To assess the displacement of the geoid caused by the atmosphere, 
we use the equation which states that it is an equipotential surface 
of gravity. Let V be the gravity potential not including the atmo- 
sphere and e6V be the internal atmospheric potential, then the geoid 
is defined by 

V(r(e), 6, $) + e6V (r(e) , 0, <fr) * const. (9) 


Here e is the ratio of the mass of the atmosphere to the mass of the 
earth plus the mass of the atmosphere, a small quantity. If e is set 
to zero we find the usual radius. Assume a Taylor expansion for r(e), 
and 

r(«) ‘ r( 0 ) ♦ (|f) E=Q * * _ 

To first order, the correction to the radius is given by 


6r = r(e) - r (0) = — 


e5V 


g * 

a result which follows from differentiation of (9). 
If we use 


M a = 4 *P 0 a H * 


then (8). becomes 


V ATM e6V 
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and hence 


6r - 


-e5V 


M 

■F r ( 


^ r ^ + T5 e 1 $ p ? ( cos 0 ) 


g m ■ id 13 -a.' 2 

To find a representative value, we set 


r - a, 


P 2 (cos 0) = 1, 


x r - m a a U + 1_ z 

5r ’ KjT a i 1 is e 


and 


The first term corresponds to a simple correction of the mass figure 
when we are operating inside the atmosphere, the second is a smaller 
latitude dependent term. To estimate their magnitude we use the 
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values on the previous pages to find 
M A 

6r 1 * jp a = S.5 meters 

X 

and 

M 

fir 2 = 5 ^ a e2) = - 5 cm - 

For present purposes, it is. apparent that a simple correction 
of the mass will suffice. The remaining latitude dependent term is 
inconsequential with radar altimeter accuracies at their present or 
anticipated levels. 

It is possible, of course, that the term 5^ above is smaller 
than the true atmospheric contribution would be if latitudinal 
density variations were taken into account. However, if the result 
is small by a factor of ten, this term will not need to be included 
in analyses until radar altimeter accuracy figures improve by a factor 
of one hundred. 
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